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ABSTRACT 


Estimates  of  the  uncertainties  attached  to  full-scale  predictions  of  submarine 
propulsion  based  on  model  tests  in  the  Large  Cavitation  Channel  (LCC)  are 
obtained  by  means  of  a  global  uncertainty  analysis.  The  analysis  takes  into  account 
all  the  component  uncertainties,  including  the  uncertainties  associated  with  the 
prediction  procedure  and  the  measurements  performed  both  at  model  scale  and  at 
full  scale,  which  influence  the  overall  uncertainty  of  full-scale  predictions. 


ADMINISTRATIVE  INFORMATION 

This  investigation  was  sponsored  by  Code  51  Program  Manager  of  the  Large 
Cavitation  Channel  Office. 


INTRODUCTION 

Estimates  of  the  imcertainties  attached  to  full-scale  predictions  of  submarine 
propulsion  based  on  propulsion  tests  in  the  LCC  are  obtained  in  this  study  by  means 
of  a  global  uncertainty  analysis.  The  analysis  takes  into  account  the  uncertainties 
associated  with  the  prediction  procedure  and  the  measurements  performed  both  at 
model  scale  and  at  full  scale.  Thus,  the  uncertainty  analysis  developed  in  the  study 
takes  into  accoimt  all  the  component  uncertainties  which  influence  the  overall 
xmcertainty  of  full-scale  predictions. 

The  prediction  procedure,  summarized  in  Appendix  A,  entails  both  tow-tank 
resistance  tests  to  determine  the  residuary-drag  coefficient,  and  propulsion  tests  in 
the  LCC.  Five  primary  model-scale  variables  are  measured  in  LCC  propulsion  tests. 
These  measured  primary  model-scale  variables  are  the  reference  velocity,  the  drag, 
and  the  propeller  rpm,  thrust  and  torque. 

The  five  measured  primary  model-scale  variables  are  used  to  determine 
several  "transformed"  model-scale  variables  by  means  of  analytical  relations.  These 
relations  are  given  in  Appendix  A.  The  transformed  model-scale  variables  include 
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four  nondimensional  variables  :  the  total-drag  coefficient,  the  advance  ratio,  the 
thrust-deduction  factor,  and  the  propulsive  efficiency. 

The  curves  representing  the  advance  ratio,  the  thrust-deduction  factor,  and 
the  propulsive  efficiency  as  functions  of  the  total-drag  coefficient  are  fundamental 
elements  of  the  model-scale  to  full-scale  extrapolation.  The  relations  used  in  this 
extrapolation  are  given  in  Appendix  A.  The  extrapolation  procedure  is  usually 
implemented  for  a  specified  full-scale  speed  or  for  a  specified  full-scale  shaft 
horsepower.  Both  cases  are  examined  in  the  global  uncertainty  analysis  developed 
in  the  study. 

The  imcertainty  analysis,  based  on  classical  expressions  for  the  errors  [1]  and 
elementary  differential  calculus,  is  expounded  in  Appendix  B.  The  Fortran-code 
implementation  of  the  expressions  for  the  uncertainties  obtained  in  Appendix  B  is 
given  in  Appendix  C.  Example  input  and  output  files  associated  with  the  Fortran- 
code  are  also  included  in  Appendix  C.  The  global  uncertainty  analysis  developed  in 
Appendices  B  and  C  provides  a  practical  tool  for  estimating  the  uncertainties  of  full- 
scale  predictions  in  terms  of  component  uncertainties  attached  to  model-scale  and 
full-scale  measurements. 

Full-scale  theoretical  predictions  are  ultimately  compared  to  values  measured 
in  full-scale  trials.  The  observed  differences  between  theoretical  predictions  obtained 
via  model-scale  tests  and  full-scale  measurements  are  usually  expressed  in  the  form 
of  a  correlation  allowance  in  the  relation  defining  the  drag  coefficient. 

The  correlation  allowance  accounts  for  aspects  of  the  full-scale  flow,  such  as 
the  hull  roughness,  that  are  not  accounted  for  in  model  tests.  The  correlation 
allowance  also  accounts  for  other  limitations  of  the  procedure  used  to  obtain  full- 
scale  predictions  from  model  tests,  notably  errors  that  are  systematically  introduced 
into  the  predictions  as  a  result  of  limitations  inherent  to  the  prediction  procedure. 
Thus,  systematic  errors  associated  with  the  characteristics  of  the  LCC  and  of  the 
experimental  set-up  used  in  the  implementation  of  the  procedure  are  largely 
included  in  the  correlation  allowance,  as  is  attested  by  the  fact  that  different 
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correlation  allowances  are  used  for  different  facilities  such  as  the  LCC  and  tow 
tanks. 

Thus,  the  correlation  allowance  largely  accounts  for  the  systematic  (bias) 
errors  associated  with  the  effects  of  the  walls  of  the  test  section  of  the  LCC,  the  strut 
holding  the  model,  the  strain  gauges,  and  the  electronic  equipment.  Therefore,  as 
long  as  no  significant  changes  are  made  in  the  characteristics  of  the  LCC,  the 
experimental  set-up  (including  the  strut,  the  strain  gauges,  and  the  electronic 
equipment)  and  the  testing  procedure,  systematic  errors  attached  to  these  aspects  of 
the  prediction  procedure  can  largely  be  ignored  in  the  uncertainty  analysis  (since 

they  are  already  included  in  the  correlation  allowance  to  a.large  extent,  as  was  noted 
previously). 

Some  errors,  however,  are  likely  to  vary  with  the  design  speed,  and  thus 
cannot  be  completely  ignored  in  the  uncertainty  analysis.  Systematic  errors  due  to 
geometrical  imprecisions  of  the  model  clearly  are  model-dependent,  and  thus 
cannot  in  principle  be  ignored  in  the  vmcertainty  analysis. 

In  summary,  it  is  appropriate  to  ignore  most  systematic  (bias)  errors  in  an 
vmcertainty  analysis  of  a  consistent  prediction  procedure  because  these  consistent 
errors  are  largely  included  in  the  correlation  allowance  attached  to  the  prediction 
procedure.  This  general  consideration  and  consideration  of  the  substantial 
difficulties  in  obtaining  realiable  estimates  of  bias  errors  —  more  precisely,  of  the 
effects  of  the  bias  errors  that  are  not  already  included  in  the  correlation  allowance  — 
®u§8®st  that  a  reasonable  practical  way  of  accounting  for  bias  errors  is  to  simply 
increase  the  precision  (random)  errors  by  means  of  a  multiplicative  factor. 
Specifically,  the  bias  errors  of  the  measured  primary  model-scale  variables  are  taken 
equal  to  the  precision  errors  of  these  variables  in  the  analysis  considered  further  on. 

The  precision  errors  attached  to  the  measured  primary  model-scale  variables 
can  be  determined  by  means  of  a  statistical  analysis  of  the  repeatability  of  model- 
scale  measurements.  This  repeatability  analysis  is  presented  in  Appendix  D. 
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Results  of  the  repeatability  analysis  presented  in  Appendix  D  and  of  the  global 
uncertainty  analysis  expounded  in  Appendix  B  are  presented  below  for  several  cases, 
with  the  purpose  of  analyzing  the  contribution  of  the  major  component 
uncertainties  which  influence  the  overall  uncertainty  of  full-scale  predictions. 

RESULTS  OF  UNCERTAINTY  ANALYSIS 

The  imcertainty  analysis  developed  in  this  study  is  applied  to  a  typical  case, 
which  is  defined  below  (and  in  the  input  file  listed  in  Appendix  C).  The  identifying 
numbers  of  the  model,  the  propeller,  and  the  resistance  (EHP)  and  propulsion  (SHP) 
tests  corresponding  to  the  case  considered  here  are 


model  no. 

propeller  no. 

EHP  exp.  no. 

SHP  exp.  no. 

XXXX 

XXXX 

XXX 

XXX 

The  water  density  and  viscosity  are 


density 

viscosity 

1.937  slug/fP 

1.084X10-^  fP/sec 

The  length  and  the  wetted  area  of  the  model,  and  the  diameter  of  the  propeller  are 


length 

area 

diameter 

22.697  ft 

138.179  fP 

0.9986  ft 

The  residuary-drag  coefficient,  determined  via  tow-tank  resistance  (EHP)  model- 
tests,  is  taken  as 

residuary-drag  coefficient 
0.00065 

The  reference  velocity,  the  propeller  rpm,  the  total  drag  Rj  ,  the  tow  force  AR,  and 
the  propeller  thrust  and  torque  in  the  propulsion  tests  are 
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ref.  velocity 

rpm 

drag 

tow  force 

thrust 

torque 

23.7  knots 

1175.0 

717.89  lbs 

191.10  lbs 

878.0  lbs 

2585.0  in-lbs 

The  slopes  of  tiie  curves  representing  the  advance  ratio,  the  thrust-deduction  factor, 
and  the  propulsive  efficiency  as  functions  of  the  total-drag  coefficient  are  equal  to 


advance  ratio 

thrust-deduction  factor 

propulsive  efficiency 

-0.249 

0.067 

-0.015 

The  length  and  the  speed  of  the  full-scale  submarine  are 


length 

speed 

380  ft 

25  knots 

Finally,  the  viscosity  of  sea  water  is  taken  as  1.282X10  ®  ft^/sec. 

As  was  already  noted,  results  of  the  repeatability  analysis  presented  in 
Appendix  D  and  the  global  uncertainty  analysis  expoimded  in  Appendices  B  and  C 
are  presented  for  several  cases  for  the  purpose  of  analyzing  the  contributions  of  the 
major  component  uncertainties  which  influence  the  overall  uncertainty  of  full- 
scale  predictions. 

Case  Ml :  contribution  of  precision-uncertainties  of  model-scale  measurements 

It  is  instructive  to  first  consider  the  full-scale  prediction-uncertainties  for  the 
case  when  only  the  imcertainties  stemming  from  propulsion  tests  in  the  LCC  are 
taken  into  accotmt.  In  this  case,  called  Ml  hereafter,  the  correlation  allowance  and 
full-scale  conditions  (i.e.  the  density  and  the  viscosity  of  sea  water,  the  geometry  of 
the  full-scale  ship  and  propeller,  the  full-scale  values  of  the  speed,  the  propeller 
rpm,  the  thrust,  the  torque,  and  the  shaft  horsepower)  are  presumed  known 
without  uncertainty.  The  residuary-drag  coefficient  (determined  via  tow-tank 
resistance  tests),  the  density  and  the  viscosity  of  the  water  in  the  LCC,  the  length  and 
the  wetted  area  of  the  model,  and  the  propeller  diameter  are  also  presumed  known 
without  imcertainty  for  case  Ml.  Furthermore,  model-scale  uncertainties  are  taken 


equal  to  the  precision  (random)  errors  determined  in  Appendix  D.  Thus,  bias  errors 
attached  to  the  measured  primary  model-scale  variables  are  not  taken  into  account 
in  case  Ml.  Case  Ml  corresponds  to  a  comparison  of  successive  model-tests  within  a 
series  of  consecutive  tests. 

Appendix  D  indicates  that  the  relative  precision  uncertainties  in  propulsion 
tests  are  approximately  equal  to  1%  for  the  reference  velocity,  0.2%  for  the  propeller 

rpm,  1.5%  for  the  drag  and  the  tow  force  AR  ,  0.5%  for  the  propeller  thrust  and 
0.3%  for  the  propeller  torque.  These  uncertainties  are  listed  below  : 

Precision  uncertainties  of  model-scale  measurements  in  propulsion  tests 


velocity 

rpm 

drag  &  tow  force 

thrust 

torque 

1% 

0.2% 

1.5% 

0.5% 

0.3% 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
previously-defined  model-scale  uncertainties  are  listed  in  the  next  table  for  two 
cases  corresponding  to  predictions  for  specified  values  of  the  full-scale  speed  or  SHP. 


Full-scale  prediction-uncertainties  for  case  Ml 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

1.0% 

2.15% 

2.6% 

2.35% 

0.0% 

SHP 

0.8% 

1.3% 

2.65% 

1.3% 

n/a 

2.35% 

The  prediction-uncertainties  listed  in  the  foregoing  table  represent  the 
contribution  of  precision  errors  of  measurements  in  the  LCC  when  all  other  sources 
of  errors  (including  bias  errors  of  model-scale  measurements,  uncertainties  of  the 
residuary-drag  coefficient  determined  via  tow-tank  resistance  tests,  uncertainties  of 
the  density  and  the  viscosity  of  the  water  in  the  LCC,  and  model-scale  geometrical 
inaccuracies)  are  ignored. 
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The  contribution  of  the  uncertainties  attached  to  the  residuary-drag 
coefficient  determined  from  tow-tank  resistance  tests,  the  density  and  the  viscosity 
of  the  water  in  the  LCC,  the  model  length  and  wetted  area,  and  the  propeller 
diameter,  are  considered  in  case  M2,  and  the  sensitivity  of  prediction-uncertainties 
attached  to  model-scale  bias  errors  is  considered  in  case  M3. 

Case  M2 :  contribution  of  uncertainties  attached  to  water  properties,  model-scale 
geometry,  and  residuary-drag  coefficient 

The  imcertainties  of  the  density  and  the  viscosity  of  the  water  in  the  LCC,  the 
model  length  and  area,  and  the  propeller  diameter  are  taken  as  is  indicated  in  the 
following  table : 


Uncertainties  of  water  properties  and  model-scale  geometry 


density 

viscosity 

length 

area 

diameter 

0.1% 

1.5% 

0.1% 

0.5% 

0.05% 

A  global  xmcertainty  analysis  of  full-scale  resistance  and  propulsion  using  tow-tank 
model  tests  [2]  shows  that  the  imcertainty  of  the  residuary-drag  coefficient  is 
approximately  8.6% ,  i.e. 

imcertainty  of  residuary-drag  coefficient 
9% 


The  prediction-imcertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
uncertainties  of  model-scale  measurements  defined  in  case  Ml  and  the 
uncertainties  of  the  residuary-drag  coefficient,  LCC-water  properties  and  model-scale 
geometry  now  considered  are  listed  in  the  next  table  for  two  cases  corresponding  to 
predictions  for  specified  values  of  the  full-scale  speed  or  SHP,  as  for  case  Ml. 
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Full-scale  prediction-uncertainties  for  case  M2 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

1.0% 

3.2% 

3.5% 

3.35% 

2.35% 

SHP 

1.15% 

1.55% 

2.8% 

1.55% 

n/a 

2.35% 

The  increase  in  uncertainties  from  case  Ml  to  case  M2  are  mainly  due  to  the 
incertainty  of  the  residuary-drag  coefficient.  In  fact,  it  can  be  verified  that  the 
uncertainties  of  the  density  and  the  viscosity  of  water  in  the  LCC,  of  the  length  and 
the  wetted  area  of  the  model,  and  of  the  diameter  of  the  propeller  are  sufficiently 
small  that  they  have  insignificant  effect  upon  the  prediction-uncertainties. 

Case  M3  :  contribution  of  model-scale  precision  and  bias  uncertainties 

As  in  cases  Ml  and  M2,  only  the  contribution  of  model-scale  uncertainties  are 
considered  in  case  M3.  Thus,  the  correlation  allowance  and  full-scale  conditions  (i.e. 
the  density  and  the  viscosity  of  sea  water,  the  geometry  of  the  full-scale  ship  and 
propeller,  the  full-scale  values  of  the  speed,  the  propeller  rpm,  the  thrust,  the 
torque,  and  the  shaft  horsepower)  are  again  presumed  known  without  uncertainty 
for  the  case  now  considered. 

As  was  already  noted  in  the  introduction,  model-scale  bias  errors  are  taken 
equal  to  the  model-scale  precision  errors  determined  in  Appendix  D  and  listed 
previously  for  cases  Ml  and  M2.  The  total  (precision  +  bias)  model-scale 
uncertainties  which  are  considered  in  case  M3  are  then  equal  to  ’2}''^  times  the 
model-scale  uncertainties  considered  in  cases  Ml  and  M2  .  Thus,  the  uncertainty  of 
the  residuary-drag  coefficient  is  now  taken  as  8.6%  X  2'''^  =  12.2%  ,  i.e. 

uncertainty  of  residuary-drag  coefficient 
12% 
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The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
previously-defined  model-scale  imcertainties  are  listed  in  the  next  table  for  two 
cases  corresponding  to  predictions  for  specified  values  of  the  full-scale  speed  or  SHP, 
as  for  cases  Ml  and  M2.  This  table  shows  that  the  prediction-imcertainties  for  case 
M3  are  approximately  equal  to  2^^^  times  the  prediction-uncertainties  for  case  M2,  as 
one  expects. 


Full-scale  prediction-imcertainties  for  case  M3 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

1.45% 

4.35% 

4.8% 

4.55% 

3.15% 

SHP 

1.55% 

2.1% 

3.9% 

2.1% 

n/a 

3.3% 

The  prediction-imcertainties  for  case  M3  may  be  regarded  as  the  uncertainties 
of  the  full-scale  predictions  obtained  using  submarine  model  testing  in  the  LCC  for 
a  specified  full-scale  submarine  and  propeller  geometry  and  specified  full-scale 
conditions.  However,  comparisons  of  full-scale  predictions  obtained  by  means  of 
model  testing  to  measurements  in  full-scale  trials  introduce  additional 
imcertainties.  These  additional  uncertainties,  called  full-scale  uncertainties 
hereafter,  stem  from  uncertainties  in  the  values  of  the  density  and  the  viscosity  of 
sea  water,  the  geometry  of  the  full-scale  ship  and  propeller,  and  the  values  of  the 
full-scale  speed,  propeller  rpm,  thrust,  torque,  and  shaft  horsepower.  The 
contribution  of  these  full-scale  uncertainties  to  the  prediction-uncertainties  is 
determined  in  case  F. 


Case  F :  contribution  of  full-scale  uncertainties 

All  model-scale  uncertainties  are  ignored  in  case  F,  which  only  considers  the 
contribution  of  full-scale  uncertainties.  Thus,  all  model-scale  variables  and  the 
correlation  allowance  are  presumed  known  without  uncertainty  in  case  F. 

The  relative  uncertainties  of  the  density  and  the  viscosity  of  sea  water,  of  the 
length  and  the  wetted-surface  area  of  the  full-scale  submarine,  and  of  the  propeller 
diameter  are  taken  as  is  indicated  in  the  following  table  : 


Uncertainties  of  full-scale  input  variables 


density 

viscosity 

length 

area 

diameter 

1% 

2% 

0.5% 

1% 

0.1% 

The  uncertainties  of  full-scale  measurements  are  considered  in  Appendix  E  .  The 
total  (precision  +  bias)  uncertainties  of  full-scale  measurements  are  taken  as 


Uncertainties  of  full-scale  measurements 


speed 

rpm 

thrust 

torque 

SHP 

0.6% 

0.4% 

3.0% 

0.9% 

0.9% 

hereafter. 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
foregoing  full-scale  uncertainties  are  listed  in  the  next  table  for  two  cases 
corresponding  to  predictions  for  specified  values  of  the  full-scale  speed  or  SHP. 


Full-scale  prediction-uncertainties  for  case  F 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

0.75% 

3.5% 

2.25% 

SHP 

0.85% 

0.7% 

3.1% 

1.2% 

0.9% 

0.9% 
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The  prediction-iincertainties  for  case  F,  which  only  considers  the  contribution 
of  full-scale  uncertainties  (with  all  other  sources  of  uncertainties  ignored),  are 
smaller  than  the  prediction-imcertainties  for  case  M3,  which  only  considers  the 
contribution  of  model-scale  uncertainties  (with  all  other  sources  of  uncertainties 
ignored). 

Case  MF :  contribution  of  model-scale  and  full-scale  uncertainties 

The  contributions  of  both  the  model-scale  uncertainties  and  the  full-scale 
vmcertainties  considered  in  cases  M3  and  F,  respectively,  are  now  combined.  The 
prediction-imcertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and  UEHP  for  the 
full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  for  this  case,  called  case  MF 
hereafter,  are  listed  in  the  next  table  for  two  cases  corresponding  to  predictions  for 
specified  values  of  the  full-scale  speed  or  SHP. 


Full-scale  prediction-uncertainties  for  case  MF 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

1.6% 

5.6% 

5.2% 

5.15% 

3.85% 

SHP 

1.75% 

2.25% 

5.0% 

2.4% 

0.9% 

3.45% 

The  uncertainties  for  case  MF  are  larger  than  the  uncertainties  for  either  case 
M3  or  case  F,  as  one  expects. 

Case  MFC  :  sensitivity  to  variations  in  the  correlation  allowance 

The  prediction-imcertainties  for  case  MF  are  based  on  the  assumption  that  the 
correlation  allowance  is  known  without  uncertainty.  However,  variations  in  the 
values  of  the  correlation  allowance  occur,  due  to  variations  in  the  full-scale 
submarine  that  are  not  accounted  for  in  model  tests  (e.g.  variations  in  the  hull 
roughness)  as  well  as  uncertainties  attached  to  both  model-scale  and  full-scale 
variables.  As  is  noted  in  the  introduction,  bias  errors  systematically  introduced  at 
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model  scale  and  full  scale  are  largely,  but  not  fully,  included  in  the  correlation 
allowance. 

The  correlation  allowance  is  taken  equal  to  0.00035  for  the  typical  case 
examined  in  the  present  uncertainty  analysis.  Experience  with  tow-tank  propulsion 
predictions  for  the  SSN  688  class  submarine  indicates  variations  of  the  correlation 
allowance  within  a  fairly  broad  range.  Inasmuch  as  the  contributions  of  model-scale 
imcertainties  and  full-scale  imcertainties  are  already  included  in  the  full-scale 
prediction-uncertainties  obtained  in  case  MF,  a  variation  of  the  correlation 
allowance  approximately  equal  to  30%  is  considered  here.  Specifically,  variations  of 
the  correlation  allowance  within  the  range 

CA  =  0.00035+7-0.0001 

are  considered  in  case  MFC.  Thus,  the  prediction-uncertainties  obtained  when  the 
effect  of  a  30%  variation  in  the  value  of  the  correlation  allowance  is  added  to  the 
model-scale  and  full-scale  imcertainties  considered  in  case  MF  is  examined  in  case 
MFC. 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  for  this  case,  called 
case  MFC,  are  listed  in  the  next  table  for  two  cases  corresponding  to  predictions  for 
specified  values  of  the  full-scale  speed  or  SHP. 


Full-scale  prediction-uncertainties  for  case  MFC 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

1.6% 

6.9% 

6.55% 

6.55% 

5.55% 

SHP 

2.25% 

2.6% 

5.15% 

2.8% 

0.9% 

3.45% 
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CONCLUSION 


In  summary,  a  tool  for  estimating  the  imcertainties  attached  to  full-scale 
predictions  of  submarine  propulsion  using  model  tests  in  the  LCC  has  been 
developed,  by  means  of  a  global  imcertainty  analysis,  and  applied  to  a  typical  case. 
The  analysis  takes  into  account  the  uncertainties  associated  with  the  prediction 
procedure  and  the  uncertainties  of  measurements  performed  both  at  model  scale 
and  at  full  scale.  Thus,  the  analysis  developed  and  applied  here  takes  into  accoimt 
all  the  component  imcertainties  which  influence  the  overall  uncertainty  of  full- 
scale  predictions. 

Estimates  of  the  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque, 
USHP  and  UEHP  attached  to  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP 
have  been  obtained  for  two  cases,  corresponding  to  predictions  for  specified  values 
of  the  full-scale  speed  or  SHP  .  Estimates  of  the  prediction-uncertainties  Uspeed, 
Urpm,  Uthrust,  Utorque,  USHP  and  UEHP  are  given  for  six  cases,  called  Ml,  M2, 

M3,  F,  MF  and  MFC. 

The  prediction-uncertainties  for  case  Ml  represent  the  contribution  of 
precision  errors  of  model-scale  measurements  in  the  LCC  when  all  other  sources  of 
errors  (including  bias  errors  of  model-scale  measurements,  uncertainty  of  the 
residuary-drag  coefficient,  uncertainties  of  the  density  and  the  viscosity  of  the  water 
in  the  LCC,  and  model-scale  geometrical  inaccuracies)  are  ignored.  Thus,  bias  errors 
attached  to  the  primary  model-scale  variables  measured  in  the  LCC  are  not  taken 
into  accoimt  in  Ml,  which  corresponds  to  successive  model  tests  within  a  series  of 
consecutive  tests. 

The  contribution  of  uncertainties  of  the  residuary-drag  coefficient,  the  density 
and  the  viscosity  of  the  water  in  the  LCC,  the  model  length  and  wetted  area,  and  the 
propeller  diameter,  are  considered  in  case  M2  .  The  increase  in  uncertainties  from 
case  Ml  to  case  M2  are  mainly  due  to  the  uncertainty  of  the  residuary-drag 
coefficient.  In  fact,  it  can  be  verified  that  the  uncertainties  of  the  density  and  the 
viscosity  of  water  in  the  LCC,  of  the  length  and  the  wetted  area  of  the  model,  and  of 
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the  diameter  of  the  propeller  are  sufficiently  small  that  they  have  insignificant 
effect. 

The  sensitivity  of  prediction-uncertainties  to  model-scale  bias  errors  is 
considered  in  case  M3.  Model-scale  bias  errors  are  taken  equal  to  the  model-scale 
precision  (random)  errors  considered  in  cases  Ml  and  M2.  Thus,  the  prediction- 
uncertainties  for  case  M3  are  equal  to  2’^^  times  the  prediction-uncertainties  for  case 
M2,  as  one  expects.  The  prediction-uncertainties  for  case  M3  may  be  regarded  as  the 
uncertainties  of  the  full-scale  predictions  obtained  using  submarine  model  testing  in 
the  LCC  for  a  specified  full-scale  submarine  and  propeller  geometry  and  specified 
full-scale  conditions. 

Comparison  of  full-scale  predictions  to  measurements  in  full-scale  trials 
introduces  additional  uncertainties.  These  additional  full-scale  uncertainties  stem 
from  uncertainties  in  the  values  of  the  density  and  the  viscosity  of  sea  water,  the 
geometry  of  the  full-scale  ship  and  propeller,  and  the  values  of  the  full-scale  speed, 
propeller  rpm,  thrust,  torque,  and  shaft  horsepower.  The  contribution  of  these  full- 
scale  uncertainties  to  the  prediction-uncertainties  is  considered  in  case  F  ,  which 
only  considers  the  contribution  of  full-scale  uncertainties  (with  all  other  sources  of 
uncertainties  ignored).  The  prediction-uncertainties  for  case  F  are  smaller  than  the 
prediction-uncertainties  for  case  M3,  which  only  considers  the  contribution  of 
model-scale  uncertainties  (with  all  other  sources  of  uncertainties  ignored). 

The  contributions  of  both  the  model-scale  uncertainties  and  the  full-scale 
uncertainties  considered  in  cases  M3  and  F  ,  respectively,  are  combined  in  case  MF  . 
Thus,  the  uncertainties  for  case  MF  are  larger  than  the  uncertainties  for  either  case 
M3  or  case  F  .  The  prediction-uncertainties  for  case  MF  are  based  on  the  assumption 
that  the  correlation  allowance  is  known  without  uncertainty. 

However,  variations  in  the  value  of  the  correlation  allowance  occur,  due  to 
variations  in  the  full-scale  submarine  that  are  not  accounted  for  in  the  model-tests 
(such  as  variations  in  the  hull  roughness),  as  well  as  uncertainties  attached  to  both 
model-scale  and  full-scale  variables.  As  is  noted  in  the  introduction,  bias  errors 
systematically  introduced  at  model  scale  and  full  scale  are  largely,  although  not 
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fully,  included  in  the  correlation  allowance.  Inasmuch  as  the  contributions  of 
model-scale  and  full-scale  imcertainties  are  already  included  in  the  full-scale 
prediction-imcertainties  evaluated  in  case  MF,  the  effect  of  a  30%  variation  in  the 
value  of  the  correlation  allowance  added  to  the  model-scale  and  full-scale 
uncertainties  considered  in  case  MF  is  examined  in  case  MFC  . 

The  cases  Ml,  M2,  M3,  F,  MF  and  MFC  are  summarized  below 


Cases  Ml,  M2,  M3,  F,  MF  and  MFC 


Ml 

only  considers  precision  imcertainties  of  model-scale  measurements 

M2 

adds  uncertainties  of  residuary-drag  coefficient,  LCC-water  properties 

and  model-scale  geometry 

M3 

considers  all  model-scale  precision  and  bias  uncertainties 

F 

only  considers  full-scale  uncertainties 

MF 

considers  all  model-scale  and  full-scale  uncertainties 

MFC 

adds  sensitivity  to  variations  in  correlation  allowance 

The  prediction-imcertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  a  specified  value  of  the  full-scale  SHP  are  listed  in  the  following  table  for 
the  six  cases  Ml,  M2,  M3,  F,  MF  and  MFC.  The  uncertainty  UPC  of  the  propulsive 
efficiency  is  also  given  in  the  table 
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Full-scale  prediction-uncertainties  for  a  specified  SHP 


case 

Uspeed 

Urpm 

Uthrust 

U  torque 

USHP 

UEHP 

UPC 

Ml 

0.8% 

1.3% 

2.65% 

1.3% 

n/a 

2.35% 

2.35% 

M2 

1.15% 

1.55% 

2.8% 

1.55% 

n/a 

2.35% 

2.35% 

M3 

1.55% 

2.1% 

3.9% 

2.1% 

n/a 

3.3% 

3.3% 

F 

0.85% 

0.7% 

3.1% 

1.2% 

0.9% 

0.9% 

n/a 

MF 

1.75% 

2.25% 

5.0% 

2.4% 

0.9% 

3.45% 

3.3% 

MFC 

2.25% 

2.6% 

5.15% 

2.8% 

0.9% 

3.45% 

3.3% 

In  summary,  it  may  be  concluded  that  the  full-scale  prediction-uncertainties 
for  a  specified  SHP  are  approximately  equal  to 

Summary  of  full-scale  prediction-uncertainties  for  a  specified  SHP 


Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

UPC 

2% 

2.5% 

5% 

3% 

1% 

3.5% 

3% 
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APPENDIX  A  ;  PREDICTION  PROCEDURE 


Primary  model-scale  variables 

Primary  model-scale  variables  are  determined  from  measurements.  Five  major  primary 
variables  are  measured  :  the  reference  velocity  V ,  the  drag  R ,  and  the  propeller  rps  n ,  thrust 
T  and  torque  Q . 


Transformed  model-scale  variables 

Transformed  model-scale  variables  are  obtained  from  the  measured  primary  variables  by 
means  of  analytical  relations.  The  major  transformed  model-scale  variables  are  the  total-drag 
coefficient  Cx ,  the  advance  ratio  Jy ,  the  thrust-deduction  factor  1  —  f  and  the  propulsive 
efficiency  Tjp .  The  relations  defining  these  transformed  variables  are  given  below. 

It  is  assumed  here  that  the  LCC  is  used  to  perform  propulsion  tests,  but  that  the  resistance 
(EHP)  tests  required  to  determine  the  residuary-drag  coefficient  Cr  are  performed  in  a  tow 
tank.  Thus,  the  uncertainty  attached  to  the  residuary-drag  coefficient  Cr  ,  determined  via  tow- 
tank  model  tests,  is  presumed  known  (i.e.,  is  an  input)  in  the  uncertainty  analysis  considered 
further  on. 

The  reference  velocity  V  for  steady  flow  past  a  model  held  fixed  inside  the  test  section  of  the 
LCC  is  determined  via  Laser  Doppler  Velocimetry  (LDV)  measurement  of  the  fluid  velocity 
at  a  section  of  the  flow  where  the  pressure  coefficient  Cp  vanishes.  The  location  of  the  Cp  =  0 
reference  section  is  determined  using  a  computational  method.  Propulsion  tests  performed  in 
the  LCC  yield  the  ideal  resistance  Ri 

Ri  =  RT-AR  (1) 


where  Rt  is  the  drag  of  the  model  without  propeller  and  AR  is  the  change  in  the  drag  due  to 
the  propeller.  The  total-drag  coefficient  Ct  of  the  model  in  a  propulsion  test  is  given  by 


Ct  — 


Ri 

pSV^/2 


(2) 


where  p  is  the  density  of  the  tank  water,  S  is  the  wetted-surface  area  of  the  model,  and  V  is 
the  previously-defined  reference  velocity.  The  advance  ratio  Jy  is  defined  as 


Jy  =  V/{nD) 


(3) 


where  n  and  D  are  the  propeller  rps  and  diameter.  The  thrust-deduction  factor  1- 1  is  given 
by 

l-t  =  Ri/T  (4) 

where  T  is  the  propeller  thrust.  The  propulsive  efficiency  rjo  is 


VRi 

2itnQ 


(5) 


where  Q  is  the  propeller  torque.  The  three  curves  representing  the  nondimensional  variables 
Jy,  1-t  and  rj£)  as  functions  of  Ct  axe  used  to  determine  full-scale  predictions. 
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Full-scale  predictions 

The  superscript  S  identifies  full-scale  variables.  No  superscript  is  used  for  model-scale 
variable.  At  a  given  ship  speed  V^,  the  total-drag  coefficient  €§• ,  the  propeller  rps  ,  thrust 
and  torque  ,  and  the  shaft  horsepower  SHP^  are  determined  using  the  relations  given 
below. 

The  friction-drag  coeflicient  Cp  is  determined  using  the  ITTC  formula 

cl  =  0.075/(C  lni?|  -  2)2  (6a) 

where  C  s;  0.4342944819 .  The  Reynolds  number  is  defined  as 

(6b) 

where  is  the  ship  length  and  is  the  kinematic  viscosity  of  sea  water.  The  total-drag 
coefficient  C|  of  the  ship  is  evaluated  using  the  relation 

C^  =  C§  +  Cr  +  Ca  (7) 


where  C|  is  the  friction-drag  coefficient,  Cr  is  the  residuary-drag  coefficient  determined  from 
resistance  (EHP)  tow-tank  model  tests,  and  the  correlation  allowance  Ca  accounts  for  differ¬ 
ences  between  the  actual  drag  coefficient  C|  and  the  predicted  drag  coefficient  C|  -|-  Cr  . 

The  propeller  rps  is  obtained  from  the  relation  » 


n 


5  = 

4DS 


(8) 


where  is  the  propeller  diameter.  Furthermore,  the  advance  ratio  Jy  is  determined  from 
the  function  Jv{Ct)  obtained  from  model  tests,  with  Ct  taken  equal  to  the  full-scale  value 
C|  predicted  by  (7) . 

The  total  drag  of  the  ship  Rj  and  the  power  required  to  overcome  R|  are  given  by 

Rf  =  C|/*“5^(K®)2/2 

EHP^  =  l4V^  =  C^  S^{V^f/2  (9) 

where  is  the  density  of  sea  water  and  is  the  wetted-surface  area  of  the  ship. 

The  thrust  exerted  by  the  propeller  is  evaluated  using  the  relation 


dS  /~>s 

s  _  -^r  _ 


psea  S^{V^f/2 


(10) 


where  the  thrust-deduction  factor  1—  is  determined  from  the  function  (1—  t)(C'T')  obtained 
from  model  tests,  with  Cp  taken  equal  to  the  full-scale  value  Cp  given  by  (7) . 

The  power  provided  to  the  propeller  is 


SHP^ 


EHP^ 


=  ^  S^{V^f/2 


(11) 
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where  the  propulsive  efficiency  rjf)  is  determined  from  the  function  obtained  from 

model  tests,  with  Ct  taken  equal  to  the  full-scale  value  given  by  (7) . 

Finally,  the  propeller  torque  is  defined  by 

=  ^  =  (12) 

27r  n*  47r  rj^ 

where  (11)  and  (8)  were  used. 

The  foregoing  relations  yield  values  of  the  shaft  horsepower  SHP^  corresponding  to  a  range 
of  values  of  the  ship  speed  .  A  plot  of  the  speed  F®  as  a  function  of  the  horsepower  SHP^ 
is  then  used  to  determine  the  ship  speed  corresponding  to  a  prescribed  value  of  the  power 
SHP^. 
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APPENDIX  B  :  GLOBAL  UNCERTAINTY  ANALYSIS 


Uncertainties  of  measured  model-scale  variables 

The  previous  relations,  which  define  model-scale  and  full-scale  variables  in  terms  of  five 
measured  primary  model-scale  variables  (reference  velocity  V ,  drag  R,  and  propeller  rps 
n ,  thrust  T  and  torque  Q) ,  can  be  used  to  determine  the  uncertainties  of  the  transformed 
model-scale  variables  and  of  the  predicted  full-scale  variables  in  terms  of  the  uncertainties 
of  the  measured  primary  variables.  These  analytical  expressions  for  the  uncertainties  of  the 
transformed  model-scale  variables  and  of  the  full-scale  predictions  are  given  below. 


Uncertainties  of  transformed  model-scale  variables 

Expression  (1)  yields 

dRi  =  dRx  ~  dAR 

The  absolute  uncertainty  of  Ri  therefore  is  given  by 


Expression  (2)  yields 

dCr  _  dRi  dp  dS  dV 

Ct  ~  Ri  p  S  ^  V 

The  relative  uncertainty  of  Ct  is  then  given  by 

Expression  (3)  defines  the  relative  uncertainty  of  the  advance  ratio  Jy  as 

fSV^ 

\Jv)~{v)'^\n)'^[D) 

The  relative  uncertainty  of  the  thrust-deduction  factor  1— t  is  defined  by  (4)  as 

The  relative  uncertainty  of  the  propulsive  efficiency  t;/)  is  defined  by  (5)  as 

In  (14)  and  (16)-(17) ,  6Ri  is  given  by  (13) . 

Uncertainties  of  full-scale  predictions 

Expressions  (6)  yield 

dC|  _  -2C  dR^  _  -2C  (dV^  dL^  du^^^\ 

Cf  ~  C\nR^-2  R^  ~  v/5:075  i/»ea  j 

-  -V^  y-yT  +  7s-  -  j 


(13) 


(14) 


(15) 


(16) 


(17) 
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where  C  =  4C^/0.075  ~  10.0593 .  We  thus  have 


dCf 

Cf 


{/«ea 


Expression  (7)  for  the  total-drag  coefficient  yields 

dCl^  =  dCf  +  dCR-{-dCA 


(18) 


which  can  be  expressed  in  the  form 


dC#  dC|  dCn  +  dCA 

nS  ~~  r^S  r^S  ^ 

We  then  have 

C#  vs) 

(19) 

where  F  and  D^''  are  defined  as 

(20) 

j^Lu_dCR  +  dCA  dLS  dv’>^'^ 

^  ^10  cl  cl 

(21) 

Expressions  (8) ,  (10) ,  (11)  and  (12)  involve  the  nondimensional  coefficients  Jv  ■,  \—t  and 
t)d  ■  These  coefficients  are  determined  from  curve  fits  (obtained  from  model  tests)  of  Jy  ,  1-t 
and  rjD  as  functions  of  Ct  ■  Let  A  stand  for  any  one  of  the  three  coefficients  Jy ,  1— t  and  rjo  ■ 
The  difference  in  the  coefficient  A  is  given  by 


dA 


cs 


dA  , 


The  first  term  represents  the  difference  in  A  at  a  given  value  of  Cj- ,  and  the  second  term 
defines  the  difference  in  A  due  to  the  uncertainty  of  .  Let  the  first  term  be  written  as  dA 
for  shortness.  The  relative  difference  in  A  may  then  be  expressed  in  the  form 


dV^\ 

j 


where  (19)  was  used.  The  relative  differences  in  the  coefficients  Jy ,  1—t  and  are  then  given 
by 


where  ct*  and  are  defined  as 


=  v/i^ 


Jv 


V  ^  \-t  dCr 


''  rjD  dCr 


(23) 


21 


Here,  expression  (20)  for  the  term  F  was  used. 


It  is  also  useful  to  define  the  notation 


dV^ 

dn^ 

^  dT^ 

r.  dQ^ 

Dp  = 

dSHP^ 

Dn  = 

J.S 

QS 

SHP^ 

(24a) 

nS  _ 

psea 

dS^ 

+  SS 

Df 

-  1-t 

II 

II 

dJy 

dD^ 

(24b) 

Expressions  (8) ,  (10) , 

(11) ,  (12) ,  (19) ,  (22)  and  (24)  yield 

Dn 

=  {1+(t'^)Dv- 

a-^D^^'-Df 

(25a) 

Dt 

=  (2- 

r  +  (T*)Dv  +  {T- 

-a*)D^''  +  D^p 

-Dt 

(25b) 

Dq 

=  (2-r 

-l-a’'- 

■<T^)Dv  +  {T-a 

^  +  a^)D^''  +  D^p-Dr,  +  Df 

(25c) 

Dp: 

=  (3- 

r-\-a'^)Dy  +  {r- 

-D, 

(25d) 

where  the  relative  differences  dJy/Jy  ,  d(l—  t^)/(l—  t^)  and  dT]f)/rif)  have  been 

taken  equal 

to  the  corresponding  model-scale  values  dJy/Jv  ,  d(l—  t)/(l—  t)  and  dtjo/flD  ■ 

The  four  relations  (25)  involve  model-scale  variables  and  differences  —  which  occur  via 
the  four  terms  dCn/CR,  dJv/Jv,  d(l  —  f)/(l~  0 1  —  and  full-scale  variables.  The 

full-scale  variables  include  the  relative  differences  ,  dD^/D^ ,  dL^/L^ 

and  dS^/S^  (which  may  be  determined  independently  and  thus  may  be  presumed  known  for 
the  purpose  of  this  uncertainty  analysis) ,  the  difference  dCA  (which  may  also  be  regarded 
as  a  given  input  for  this  analysis) ,  and  the  five  terms  dV^/V^y  dn^/n^,  dT^/T^,  dQ^ /Q^, 
dSHP^ /SHP^ .  Thus,  four  of  these  five  terms  may  be  determined  from  any  one  of  them.  Specif¬ 
ically,  the  four  relations  (25) ,  which  define  the  relative  differences  dn^/n^,  dT^/T^,  dQ^ /Q^ 
and  dSHP^/SHP^  in  terms  of  dV^/V^,  can  be  expressed  in  four  alternative  forms  which  define 
the  full-scale  prediction  uncertainties  in  terms  of  dn^/n®,  dT^/T^,  dQ^/Q^  or  dSHP^/SHP^. 
These  four  alternative  forms  are  considered  in  [2]  for  the  similar  uncertainty  analysis  of  full- 
scale  submarine  propulsion  predictions  using  tow-tank  model  tests.  Only  the  most  useful 
alternative  form  of  the  relations  (25) ,  which  defines  the  relative  differences  dV^/V^,  dn^ /n^, 
dT^/T^  and  dQ^/Q^  in  terms  of  dSHP^/SHP^ ,  is  considered  here. 

This  alternative  form  of  the  relations  (25)  is 

{3  -  T  +  a^)  Dv  =  Dp  -  {r-  +  Dr,  (26a) 

(3-T  +  a'^)DN  =  il+(T''){Dp-D^^  +  Dr,)-{T-(T^  +  3a'’)D^'' 

-{3-T  +  cr^)Df  (26b) 

(Z-T  +  a^)DT  =  i2-T  +  a^){Dp  +  D„)  +  {r+2<T'^ -3a^)D^'' 

+  {l+a^ -a*)D^-{3-r  +  (T'^)Dt  (26c) 

{3-T  +  a^)DQ  =  {2-r  +  (T^ -(T-’)Dp  +  {r-(T'^  +  3a-^)D^'' 

+  {l+a-^){Df-Dr,)  +  {3-r  +  <T^)D^  (26d) 
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As  was  already  noted,  the  foregoing  relations  define  the  relative  differences  dV^/V^,  dn^/n^, 
dT^/T^  and  dQ^ /Q^  in  terms  of  dSHP^/SHP^  . 

Expressions  (9)  and  (11)  show  that  the  relative  diflference 


„  _dEHP^ 
EHPS 


(27) 


can  be  obtained  from  the  foregoing  expressions  for  the  relative  difference  Dp  =  dSHP^/SHP^ . 
Specifically,  expression  (25d)  3delds 

DE  =  {Z-r)Dv  +  rD^‘'  +  D^p  (28a) 

The  relations  EHP^  =  r]D  SHP^ ,  (22)  and  (26a)  yield 

iS-r  +  (7^)DE  =  {3-T)iDp  +  Dr,)  +  a^3D^''  +  D^p)  (28b) 


The  relative  uncertainties  6V^/V^,  6n^/n^,  6T^/T^,  SQ^/Q^  and  SSHP^/SHP^  corre¬ 
sponding  to  the  relative  diflFerences  dV^/V^,  dn^/n^,  dT^/T^,  dQ^ /Q^  and  dSHP^/SHP^ 
defined  in  the  foregoing  alternative  relations  are  readily  determined  by  taking  the  square  root 
of  the  sum  of  the  square  of  every  term  in  these  relations.  Thus,  we  define  the  notation 


Vv 


10(C'|)3  j  ■’'I  j 

■(SJ  "■ 


Up 


“  \  / 


+ 


(?) 


Up  = 


Un 


fSSHP^^ 


=  gj 


Uq 


Ue 


fSEHP^^ 
^EHPS ) 


Ut 


SHPS  ) 


■(¥)’ 


-ml 


(29a) 

(29b) 

(29c) 

(29d) 


corresponding  to  (21)  and  (24) .  We  also  define  the  relative  uncertainties  attached  to  the 
full-scale  measurements  of  V^,  n®,  T^,  ,  SHP^  and  EHP^,  i.e. 


where  the  subscript  or  superscript  fsm  means  full-scale  measurement.  The  uncertainty 
attached  to  the  effective  horsepower  EHP^  of  the  ship  is  not  defined  in  practice  because 
measurements  of  EHPff^  are  not  available.  Thus,  the  term  EHP^f^  may  be  ignored  in  the 
expressions  given  below. 
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Expressions  (25)  and  (28a)  yield 

UN  =  il+(r-^f  )2  U^''  +  f/f  +  t/j?’”  (30a) 

[/t  =  (2  -  r  +  )2 +  (r-  (7‘  )2 uS^  +  +  Ut  +  (30b) 

C/q  =  ( 2  -  r  +  -  (7-^  )2  +  ( r-  a"  +  (T^  f  Uh''  +  +  Ur,  +  Uf  +  (30c) 

t/p  =  ( 3  -  r  +  a"  )2 +  ( r-  a” )2  uh''  +  +  Ur,  +  U^p'^  (30d) 

=  ( 3  -  r  )2  +  r2  uh''  +  [/p  +  (30e) 

Similarly,  (26)  and  (28b)  yield 

( 3  -  r  +  a’'  )2  C7v  =  uh'”"  +  {T  -  Uh''  +  U^  + Ur,  +  {3 -T  +  cr^f  Ul^^  (31a) 

{3-T  +  (T^fUN  =  {l+<T^f{Ufp"^  +  U^  +  Ur,)  +  iT-a^  +  3a'^)^Uh'' 

+  {3-r  +  (T^)^Uf  +  Ug'^)  (31b) 

{3-T  +  (r^fUT  =  {2-r  +  (T*  )2(f7/®"’  +  [/„)  +  (r+2<T’'-3CT‘  fuh'' 

+  {l+<T^-<T^)^U^  +  {3-r  +  a'^)^{Ut  +  Uf'^)  (31c) 


( 3  -  r  +  a’' )2  t/g  =  ( 2  -  r  +  (t’' -  (T-' )2  c//**"  +  ( r- (r'' +  3  )2 

+  (l+<T-^)2(t/p«  +  t/r,)  +  (3-r  +  (r’')2(f/f +  c;^®"’)  (31d) 

{3-r  +  cT'^)^UE  =  {3-T)'^{Up'^  +  Ur,)  +  {a'^)H9Uh‘'  +  U^) 

+  i3-T  +  <T^fuh""'  (31e) 

The  two  sets  of  alternative  expressions  (30)  and  (31) ,  and  expressions  (29) ,  (23)  and  (20) , 
define  the  uncertainties  attached  to  the  full-scale  predictions  of  the  ship  speed  V^,  the  propeller 
rpm  N^,  thrust  and  torque  Q^,  the  shaft  horsepower  SHP^  and  the  effective  horsepower 
EHP^  for  the  two  cases  in  which  or  SHP^  are  held  constant  (within  the  accuracy  of 
full-scale  measurements) . 
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APPENDIX  C : 

FORTRAN-CODE  &  INPUT-OUTPUT  FILES 

The  source  file  LCCUAi  of  the  program  LCCUA  ,  which  represents  the 
Fortran  implementation  of  the  imcertainty  analysis  expounded  in  Appendix  B  ,  is 
given  in  Appendix  C.  The  symbols  defined  in  the  analysis  and  used  in  the  Fortran- 
code  LCCUA  are  fairly  consistent. 

An  example  of  the  input  file  LCCUA.in  required  by  LCCUA.f ,  and  of  the 
corresponding  output  file  LCCUA.out  generated  by  LCCUA.f ,  is  also  given  in  this 
Appendix.  The  attached  example  input  file  LCCUA.in  and  output  file  LCCUA.out 
corresponds  to  the  previously-defined  case  MFC  ,  in  which  the  uncertainties 
attached  to  model-scale  and  full-scale  variables  and  to  the  value  of  the  correlation 
coefficient  are  included. 
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FORTRAN  CODE 


c  Global  uncertainty  analysis  of  full-scale  submarine 
c  propulsion  predictions  using  model  tests  in  LCC 
c  Francis  Noblesse  (April  98) 
c 

program  LCCUA 
c 

character  LCCxp*50  ,  date*50  ,  model*50  ,  prop*50  , 

&  EHPxp*50  ,  SHPxp*50  ,  comment*80 
c 

real  ro ,  nu  ,  Uro ,  Unu ,  L ,  S ,  D ,  UL ,  US ,  UD , 

&  CR ,  UCR ,  Vknot ,  RT ,  DelR ,  nrpm  ,  T ,  Qinlb  , 

&  UV ,  URT ,  UDelR ,  Un ,  UT ,  UQ ,  JVCT ,  tdCT ,  etaCT , 
&  nusea  ,  Uros  ,  Unus  ,  LS  ,  VSknot ,  ULS  ,  USS  ,  UDS  , 

&  UVfsm  ,  UNfsm  ,  UTfsm  ,  UQfsm  ,  UPfsm  ,  CA  ,  dCA  , 
&  V ,  VS ,  n ,  Q ,  CTCR ,  CFCR ,  Ri ,  CT ,  CF , 

&  JV ,  td  ,  eta  ,  Uro2  ,  Unu2  ,  UL2  ,  US2 ,  UD2 , 

&  UV2  ,  URi ,  URi2 ,  Un2  ,  UT2  ,  UQ2  ,  UCA , 

&  UCTCR2  ,  UCTCR ,  UCFCR2  ,  UCFCR , 

&  UCT2 ,  UCT ,  UCF2 ,  UCF , 

&  UJV2  ,  Utd2  ,  Ueta2,  UJV  ,  Utd  ,  Ueta  , 

&  CFS  ,  CTS ,  UroS2  ,  UJD2  ,  cofCFS ,  UCLnu2  , 

&  Gamma  ,  sigmaj ,  sigmat ,  sigmeta  , 

&  UVfsm2  ,  UNfsm2  ,  UTfsm2  ,  UQfsm2  ,  UPfsm2  , 

&  UN2V  ,  UT2V  ,  UQ2V  ,  UP2V  ,  UE2V  , 

&  UW  ,  UNV  ,  UTV  ,  UQV  ,  UPV  ,  UEV  ,  UAV  , 

&  UV2P  ,  UN2P ,  UT2P ,  UQ2P ,  UE2P , 

&  UPP  ,  UVP  ,  UNP  ,  UTP ,  UQP ,  UEP ,  UAP 
c 

c  READ  INPUT  VARIABLES 
c 

open(ll,file='LCCUA.in’,status='old') 

c 

read(ll,*)  LCCxp 
read(ll/)  date 
read(ll/)  model 
read(ll/)  prop 
read(ll/)  EHPxp 
read(ll/)  SHPxp 
read(ll/)  comment 
read(ll,*) 
read(ll/) 
read(ll,*) 
read(ll,*) 
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read(ll/) 
read(ll,*)  ro ,  nu 
read(ll,*) 
read(ll/) 
read(ll,*) 

read(ll/)  Uro,  Unu 

read(ll/) 

read(ll/) 

read(ll,*) 

read(ll/)L,S,D 

read(ll,*) 

read(ll/) 

read(ll,*) 

read(ll/)  UL ,  US ,  UD 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll/)  CR 

read(ll,*) 

read(ll,*) 

read(ll/)  UCR 

read(ll/) 

read(ll/) 

read(ll,*) 

read(ll,*) 

read(ll/)  Vknot ,  RT  ,  DelR 

read(ll/) 

read(ll,*) 

read(ll/)  nrpm ,  T ,  Qinlb 

read(ll,*) 

read(ll/) 

read(ll,*) 

read(ll/)  UV  ,  URT  ,  UDelR ,  Un  ,  UT  ,  UQ 

read(ll,*) 

read(ll,*) 

read(ll/) 

read(ll/)  JVCT ,  tdCT ,  etaCT 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll/) 

read(ll,*) 

read(ll/)  nusea 

read(ll,*) 

read(ll,*) 
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read(ll,*) 

read(ll,*)  Uros  ,  Unus 
read(ll,*) 
read(ll/) 
read(ll,*) 

read(ll/)  LS ,  VSknot 
read(ll/) 
read(ll,*) 
read(ll,*) 

read(ll/)  ULS  ,  USS  ,  UDS 

read(ll/) 

read(ll/) 

read(ll,*) 

read(ll/)  UVfsm  ,  UNfsm  ,  UTfsm  ,  UQfsm  ,  UPfsm 

read(ll/) 

read(ll/) 

read(ll/) 

read(ll,*)  CA ,  dCA 
c 

close(ll,status=’keep') 

c 

c  About  notation  : 

c  U  stands  for  relative  Uncertainty 

c  fsm  stands  for  Full-Scale  Measurement  uncertainty 

c 

c  PRELIMINARY  TRANSFORMATIONS 

c 

c  Rescale  nu  and  nusea 
c 

nu  =  nu  /  100000. 
nusea  =  nusea  /  100000. 
c 

c  Transform  speeds  from  knots  to  ft/sec 

c 

V  =  1.6878  Vknot 
VS  =  1.6878  *  VSknot 
c 

c  Transform  rpm  into  rps 

c 

n  =  nrpm  /  60. 
c 

c  Transform  torque  from  in-lb  to  ft-lb 

c 

Q  =  Qinlb/12. 
c 

c  Transform  input  percent  uncertainties 
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c 

Uro  =  0.01  *  Uro 
Unu  =  0.01  *  Unu 
c 

UL  =  0.01  *  UL 
US  =  0.01  *  US 
UD  =  0.01  *  UD 
c 

UCR  =  0.01  *  UCR 
c 

UV  =  0.01  *  uv 
URT  =  0.01  *  URT 
UDelR  =  0.01  *  UDelR 
c 

Un  =  0.01  *  Un 
UT  =  0.01  *  UT 
UQ  =  0.01  UQ 
c 

Uros  =  0.01  *  Uros 
Unus  =  0.01  *  Unus 
c 

ULS  =  0.01  ULS 
USS  =  0.01  *  USS 
UDS  =  0.01  *  UDS 
c 

UVfsm  =  0.01  *  UVfsm 
UNfsm  =  0.01  *  UNfsm 
UTfsm  =  0.01  *  UTfsm 
UQfsm  =  0.01  UQfsm 
UPfsm  =  0.01  UPfsm 
c 

c  MODEL-SCALE  VARIABLES 

c 

c  Compute  Ri ,  CT  &  CF 
c 

Ri  =  RT  -  DelR 

CT  =  2.  *  Ri  /  (  ro  *  S  *  V V  ) 

CF  =  0.075  /  ( LOG10(  V  *  L  /  nu )  -  2.  )**2 
c 

c  Compute  JV  ,  td=l-t  &  eta=etaD 

c 

JV  =  V/(n*D) 
td  =  Ri  /  T 

eta  =  V  *  Ri  /  ( 6.2831853  *  n  Q ) 
c 

c  PRELIMINARY  CALCULATIONS  FOR  MODEL-SCALE  UNCERTAINTIES 
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c 


c 


c 

c 


c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


c 

c 

c 

c 

c 


Uro2  =  Uro  Uro 
Unu2  =  Unu  *  Unu 

UL2  =  UL  *  UL 
US2  =  US  *  US 
UD2  =  UD  *  UD 

UV2  =  UV  UV 

URi  =  (  URT  *  RT  )**2  +  (  UDelR  *  DelR  )**2 
URi  =  SQRT(  URi )  /  Ri 
URi2  =  URi  *  URi 


Un2  =  Un  *  Un 
UT2  =  UT  *  UT 
UQ2  =  UQ  *  UQ 

MODEL-SCALE  UNCERTAINTIES 
Compute  dCT  /  CT 

UCT2  =  URi2  +  Uro2  +  US2  +  4.  *  UV2 
UCT  =  100.  *  SQRT(  UCT2  ) 


Compute  dCF  /  CF 

UCF2  =  10.  *  CF  *  (  UV2  +  UL2  +  Unu2  ) 
UCF  =  100.  *  SQRT(  UCF2  ) 


Compute  dJV  /  JV ,  dtd  /  td  &  deta  /  eta 

UJV2  =  UV2  +  Un2  +  UD2 

Utd2  =  URi2  +  UT2 

Ueta2  =  UV2+  URi2  +  Un2  +  UQ2 

UJV  =  100.  *  SQRT(  UJV2  ) 

Utd  =  100.  *  SQRT(  Utd2  ) 

Ueta  =  100.  *  SQRT(  Ueta2  ) 

PRELIMINARY  CALCULATIONS  FOR  FULL-SCALE  UNCERTAINTIES 
Compute  CFS  &  CTS 

CFS  =  0.075  /  ( LOG10(  VS  *  LS  /  nusea  )  -  2.  )**2 
CTS  =  CFS  +  CR  +  CA 
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c 

c 

c 


c 

c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


c 


c 


Compute  (drhosea/rhosea)'^2+(dSS/SS)^2  &  (dJV/JV)^2+(dDS/DS)^2 

UroS2  =  Uros'^*2  +  USS**2 
UJD2  =  UJV2  +  UDS**2 

Compute  UCLnu2  &  Gamma 

cofCFS  =  10.  *  CFS*’'3 

UCLnu2  =  (  CR  *  UCR  )**2  +  dCA**2 
UCLnu2  =  UCLnu2  /  cofCFS  +  ULS**2  +  Unus**2 

cofCFS  =  SQRT(  cofCFS ) 

Gamma  =  cofCFS  /  CIS 

Compute  sigmaj ,  sigmat  &  sigmeta 

sigmaj  =  cofCFS  *  JVCT  /  JV 
sigmat  =  cofCFS  *  tdCT  /  td 
sigmeta  =  cofCFS  *  etaCT  /  eta 


Compute  squares  of  full-scale  measurement  uncertainties 

UVfsm2  =  UVfsm  *  UVfsm 
UNfsm2  =  UNfsm  *  UNfsm 
UTfsm2  =  UTfsm  UTfsm 
UQfsm2  =  UQfsm  *  UQfsm 
UPfsm2  =  UPfsm  *  UPfsm 

FULL-SCALE  UNCERTAINTffiS  @  a  GIVEN  SPEED 

UN2V  =  UVfsm2  *  (  1.  +  sigmaj  )**2  +  UCLnu2  *  sigmaj’**2 
UN2V  =  UN2V  +  U7D2  +  UNfsm2  ^ 

UT2V  =  UVfsm2  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2V  =  UT2V  +  UCLnu2  (  Gamma  -  sigmat  )**2 
UT2V  =  UT2V  +  UroS2  +  Utd2  +  UTfsm2 


UQ2V  =  UVfsm2  *  (2.-  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2V  =  UQ2V  +  UCLnu2  *  (  Gamma  -  sigmeta  +  sigmaj  )**2 
UQ2V  =  UQ2V  +  UroS2  +  Ueta2  +  UJD2  +  UQfsm2  ° 


UP2V  =  UVfsm2  *  (  3.  -  Gamma  +  sigmeta  )**2 
UP2V  =  UP2V  +  UCLnu2  (  Gamma  -  sigmeta  )’''’'’2 
UP2V  =  UP2V  +  UroS2  +  Ueta2  +  UPfsm2 
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UE2V  =  UVfsm2  *  (  3.  -  Gamma  )**2  +  UCLnu2  *  Gamma**2  +  UroS2 
c 

UW  =  100.  *  UVfsm 
UNV  =  100.  SQRT(  UN2V  ) 

UTV  =  100.  *  SQRT(  UT2V  ) 

UQV  =  100.  *  SQRT(  UQ2V  ) 

UPV  =  100.  *  SQRT(  UP2V  ) 

UEV  =  100.  *  SQRT(  IJE2V  ) 

UAV  =  20.  *  SQRT(  UVfsm2  +  UN2V  +  UT2V  +  UQ2V  +  UP2V  ) 
c 

c  FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  SHAFT  HORSEPOWER 
c 

UV2P  =  UPfsm2  +  UCLnu2  *  (Gamma-sigmeta)**2  +  UroS2  +  Ueta2 
UV2P  =  UV2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  UVfsm2 
c 

UN2P  =  (  LJPfsm2  +  UroS2  +  Ueta2  )*(!.  +  sigmaj  )**1 
UN2P  =  LJN2P  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**1 
UN2P  =  UN2P  /  (  3.  -  Gamma  +  sigmeta  )**1  +  UJD2  +  UNfsm2 
c 

UT2P  =  (  UPfsm2  +  Ueta2  )  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2P  =  UT2P  +  UCLnu2  *  (Gamma+2.*sigmeta-3.*sigmat)**2 
Lnr2P  =  UT2P  +  UroS2  *  ( 1.  +  sigmeta  -  sigmat  Y*2 
UT2P  =  UT2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  Utd2  +  UTfsm2 
c 

UQ2P  =  LJPfsm2  *  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2P  =  UQ2P  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  sigmaj  )**2 
UQ2P  =  UQ2P  +  (  UroS2  +  Ueta2  )  M  1-  +  sigmaj  y*2 
UQ2P  =  UQ2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  UJD2  +  UQfsm2 
c 

UE2P  =  (  UPfsm2  +  Ueta2  )  *  (  3.  -  Gamma  )**2 
UE2P  =  UE2P  +  (  9.  *  UCLnui  +  UroS2  )  *  sigmeta*’^2 
UE2P  =  UE2P  /  (  3.  -  Gamma  +  sigmeta  )**2 
c 

UPP  =  100.  *  UPfsm 
UVP  =  100.  *  SQRT(  UV2P  ) 

UNP  =  100.  *  SQRT(  UN2P  ) 

UTP  =  100.  *  SQRT(  UT2P  ) 

UQP  =  100.  SQRT(  UQ2P ) 

UEP  =  100.  *  SQRT(  UE2P  ) 

UAP  =  20.  *  SQRT(  UV2P  +  UN2P  +  UT2P  +  UQ2P  +  UPfsm2  ) 
c 

c  WRITE  INPUT  VARIABLES  &  OUTPUT  RESULTS 
c 

c  Express  relative  uncertainties  in  percent 
c 


32 


Uro  =  100.  *  Uro 
Unu  =  100.  *  Unu 
c 

UL  =  100.  *  UL 
US  =  100.  US 
UD  =  100.  *  UD 
c 

UCR  =  100.  *  UCR 
c 

UV  =  100.  *  UV 
URT  =  100.  *  URT 
UDelR  =  100.  *  UDelR 
c 

Un  =  100.  *  Un 
UT  =  100.  *  UT 
UQ  =  100.  *  UQ 
c 

Uros  =  100.  *  Uros 
Unus  =  100.  *  Unus 
c 

ULS  =  100.  *  ULS 
USS  =  100.  *  USS 
UDS  =  100.  *  UDS 
c 

UVfsm  =  100.  *  UVfsm 
UNfsm  =  100.  *  UNfsm 
UTfsm  =  100.  *  UTfsm 
UQfsm  =  100.  *  UQfsm 
UPfsm  =  100.  *  UPfsm 
c 

UCA  =  100.  *  dCA  /  CA 
c 

open(12,file='LCCUA.out',status='new') 

c 

write(12/)  LCCxp 
write(12,*)  date 
write(12/)  model 
write(12,*)  prop 
write(12/)  EHPxp 
write(12,*)  SHPxp 
write^2,*)  comment 
c 

write(12/) 

write(12,*)  '  INPUT  VARIABLES’ 
write(12,*) 
c 


33 


write(12,*) 

write(12/)  ’  TANK-WATER  PROPERTIES’ 

write(12,*) 

write(12,101)  ro 

write(12,102)  nu 

write(12,103)  Uro 

write(  12,104)  Unu 

write(12,’') 

write(12,*)  '  MODEL  GEOMETRY’ 

write(12,’^) 

write(12,105)  L 

write(12,106)  S 

write(12,107)  D 

write(12,108)  UL 

write(12,109)  US 

write(12,110)  UD 

write(12,’^) 

write(12,*)  ’  MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS’ 

write(12,*) 

write(12,lll)  CR 

write(12,112)  UCR 

write(12,*) 

writeCU,-^)  ’  MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS' 

write(12,*) 

write(12,114)  Vknot 

write(12,115)  RT 

write(12,116)  DelR 

write(12,117)  nrpm 

write(12,118)  T 

write(12,119)  Qinlb 

write(12,’*) 

write(12,*)  ’  UNCERTAINTIES  OF  MODEL-SCALE  MEASUREMENTS’ 

write(12,’^) 

write(12,120)  UV 

write(12,121)  URT 

write(12,122)  UDelR 

write(12,123)  Un 

write(12,124)  UT 

write(12,125)  UQ 

write(12,’^) 

write(12,’^)  ’  OTHER  MODEL-SCALE  VARIABLES’ 
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write(12/) 
write(12,126)  JVCT 
write(12,127)  tdCT 
write(12,128)  etaCT 

write(12/) 

write(12/)  '  SEA-WATER  PROPERTIES' 

write(12,*) 

write(12,151)  nusea 

write(12,152)  Uros 

write(12,153)  Unus 

write(12,*) 

write(12/)  ’  FULL-SCALE  SHIP' 
write(12/) 
write(12454)  LS 
write(12,155)  VSknot 
write(12,156)  ULS 
write(12457)  USS 
write(12458)  UDS 

write(12,*) 

write(124)  '  UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS' 

write(124) 

write(12,159)  UVfsm 

write(12,160)  UNfsm 

write(12,161)  UTfsm 

write(12,162)  UQfsm 

write(12,163)  UPfsm 

write(124) 

write(124)  '  SCALING  ALLOWANCE’ 
write(124) 
write(12,171)  CA 
write(12,172)  dCA 

write(12,*) 

write(124) 

write(124)  '  OUTPUT  VARIABLES’ 
write(124) 

write(124) 

write(124)  ’  MODEL-SCALE  VARIABLES’ 

write(124) 

write(12,204)  CT 

write(12,205)  CF 
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write(12,206)  CR 
write(12,*) 
write(12,207)  JV 
write(12,208)  td 
write(12,209)  eta 

write(12/) 

write(12/)  ’  UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES' 

write(12/) 

write(12,213)  UCT 

write(12,214)  UCF 

write(12,215)  UCR 

write(12,216)  UCA 

write(12,*) 

write(12,217)  UJV 

write(12,218)  Utd 

write(12,219)  Ueta 

write(12,*) 

write(12/)  '  FULL-SCALE  VARIABLES' 

write(12/) 

write(12,220)  CTS 

write(12,221)  CFS 

write(12,222)  CR 

write(12,223)  CA 

write(12,*) 

write(12/) '  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 

SPEED’ 

write(12,*) 

write(12,311)  UW 

write(12,312)  UNV 

write(12,313)  UTV 

write(12,314)  UQV 

write(12,315)  UPV 

write(12,316)  UEV 

write(12,317)  UAV 

write  (12,*) 

write(12,*) '  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 
SHP’ 

write(12,*) 
write(12,321)  UVP 
write(12,322)  UNP 
write(12,323)  UTP 
write(12,324)  UQP 
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write(12,325)  UPP 
write(12326)  UEP 
write(12327)  UAP 
c 

close(12,status='keep') 

c 

c  FORMATS 
c 

101  formate  water  density  (slug/ft**3) :  ',F11.3) 

102  formate  kinematic  viscosity  coefficient  (ft**2/sec)  :  ',E11.4) 

103  formate  percent  uncertainty  of  density  :  ’,F8.2) 

104  formate  percent  xmcertainty  of  viscosity  :  ',F8.2) 

c 

105  formate  length  of  model  (ft) :  ’,F8.3) 

106  formate  wetted  area  (ft**2) :  ',F8.3) 

107  formate  diameter  of  propeller  (ft) :  ’/F8.4) 

108  formate  percent  imcertainty  of  length  :  ',F8.2) 

109  formate  percent  imcertainty  of  wetted  surface  :  ',F8.2) 

110  formate  percent  uncertainty  of  prop  diameter  :  ',F8.2) 

c 

111  formate  residuary-resistance  coefficient :  ',F10.5) 

112  formate  percent  uncertainty  of  residuary-resistance  coef :  ',F8.2) 
c 

114  formate  reference  velocity  (knots)  :  ',F8.2) 

115  formate  drag  RT  (lbs) :  ',F8.2) 

116  formate  drag  DeltaR  (lbs) :  ’,F8.2) 

117  formate  propeller  rpm :  ',F8.2) 

118  formate  propeller  thrust  (lbs) :  ',F8.2) 

119  formate  propeller  torque  (in-lbs)  :  ’,F8.2) 
c 

120  formate  percent  uncertainty  of  reference  velocity  (knots)  :  ',F8.2) 

121  formate  percent  uncertainty  of  drag  RT  (lbs) :  ',F8.2) 

122  formate  percent  uncertainty  of  drag  DeltaR  (lbs) :  ’,F8.2) 

123  formate  percent  imcertainty  of  propeller  rpm  :  ',F8.2) 

124  formate  percent  uncertainty  of  prop  thrust  (lbs) :  ',F8.2) 

125  formate  percent  uncertainty  of  prop  torque  (in-lbs) :  ’,F8.2) 

c 

126  formate  slope  d  JV  /  d  CT  :  ’,F9.3) 

127  formate  slope  d  (1-t)  /  d  CT :  ’,F9.3) 

128  format('  slope  d  etaD  /  d  CT  :  ’,F9.3) 
c 

151  formate  kinematic  viscosity  coefficient  (ft**2/sec)  :  ',E11.4) 

152  formate  percent  uncertainty  of  density  :  ',F8.2) 

153  formate  percent  uncertainty  of  viscosity  :  ',F8.2) 
c 

154  formate  ship  length  (ft) :  ’,F8.2) 
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155  formate  ship  speed  (knots) :  ',F11.1) 

c 

156  formate  percent  uncertainty  of  ship  length  :  ’,F8.2) 

157  formate  percent  uncertainty  of  wetted  surface  :  ',F8.2) 

158  formate  percent  uncertainty  of  prop  diameter  :  ',F8.2) 

c 

159  formate  percent  uncertainty  of  ship  speed  :  ',F8.2) 

160  formate  percent  uncertainty  of  prop  rpm  :  ',F8.2) 

161  formate  percent  uncertainty  of  prop  thrust :  ',F8.2) 

162  formate  percent  uncertainty  of  prop  torque  :  ',F8.2) 

163  formate  percent  uncertainty  of  SHP  :  ',F8,2) 

c 

171  formate  correlation  allowance  :  ',F11.5) 

172  formate  uncertainty  of  allowance  :  ',F11.5) 
c 

201  formate  total  resistance  coefficient  CT  :  ’,F10.5) 

202  formate  friction  resistance  coefficient  CF  :  ',F10.5) 

203  formate  residuary  resistance  coefficient  CR  :  ',F10.5) 

c 

204  formate  total  resistance  coefficient  CT  :  ',F10.5) 

205  formate  friction  resistance  coefficient  CF  ;  ’,F10.5) 

206  formate  residuary  resistance  coefficient  CR  :  ’,F10.5) 

c 

207  formate  advance  ratio  JV  :  ',F8.2) 

208  formate  thrust-deduction  factor  1-t :  ',F8.2) 

209  formate  propulsive  efficiency  etaD  :  ',F8.2) 
c 

210  formate  percent  uncertainty  of  CT  :  ',F8.2) 

211  formate  percent  uncertainty  of  CF  :  ’,F8.2) 

212  formate  percent  uncertainty  of  CR  :  ’,F8.2) 

c 

213  formate  percent  uncertainty  of  CT  :  ',F8.2) 

214  formate  percent  uncertainty  of  CF  :  ',F8.2) 

215  formate  percent  uncertainty  of  CR  :  ',F8.2) 

216  formate  percent  uncertainty  of  CA  :  ',F8.2) 

c 

217  formate  percent  uncertainty  of  JV  :  ',F8.2) 

218  formate  percent  uncertainty  of  1-t :  ’,F8.2) 

219  formate  percent  uncertainty  of  etaD  :  ’,F8.2) 

c 

220  formate  total  resistance  coefficient  CT  :  ',F10.5) 

221  formate  friction  resistance  coefficient  CF  :  ’,F10.5) 

222  formate  residuary  resistance  coefficient  CR  :  ',F10.5) 

223  formate  correlation  allowance  coefficient  CA  :  ',F10.5) 
c 

311  formate  percent  uncertainty  of  ship  speed  :  ',F8.2) 
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312  fonnatC  percent  uncertainty  of  prop  rpm  :  ',F8.2) 

313  fonnat('  percent  uncertainty  of  prop  tiirust :  ',F8.2) 

314  format('  percent  uncertainty  of  prop  torque  :  ',F8.2) 

315  formate  percent  uncertainty  of  SHP  :  ',F8.2) 

316  formate  percent  uncertainty  of  EHP :  ',F8.2) 

317  formate  percent  overall  imcertainty  :  ',F8.2) 

c 

321  formate  percent  imcertainty  of  ship  speed  :  ',F8.2) 

322  formate  percent  imcertainty  of  prop  rpm  :  ',F8.2) 

323  formate  percent  uncertainty  of  prop  thrust :  ’,F8.2) 

324  formate  percent  uncertainty  of  prop  torque  :  '/F8.2) 

325  formate  percent  uncertainty  of  SHP  :  'J’8.2) 

326  formate  percent  uncertainty  of  EHP  :  ',F8.2) 

327  formate  percent  overall  uncertainty  :  ',F8.2) 

c 

stop 

end 
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EXAMPLE  INPUT  FILE 


’  LCC  EXP ' 

'  DATE:  XXXX’ 

'  MODEL  No.  XXXX’ 

’  PROPELLER  No.  XXXX’ 

’  EHP  EXPERIMENT  No.  XXX  ' 

’  SHP  EXPERIMENT  No.  XXX  ’ 

’  COMMENTS:  CASE  MFC  ’ 

MODEL-SCALE  VARIABLES  AND  UNCERTAINTIES 

Tank-water  density  (rho)  and  kinematic  viscosity  (nu) 
rho  (slug/ft**3)  nu  X  10*’^5  (ft**2/sec) 

1.9367,  1.084 

Percent  relative  uncertainties  of  rho  &  nu 
density  kinematic  viscosity 
0.14 ,  2.1 

Model  geometry 

length  (ft)  area  (ft**2)  prop  diameter  (ft) 

22.697,  138.179  ,  0.9986 

Percent  relative  uncertainties  of  model  geometry 
length  area  prop  diameter 
0.14 ,  0.71 ,  0.07 

RESISTANCE  (EHP)  TESTS 

Residuary-resistance  coefficient 
0.00065 

Percent  relative  uncertainty  of  residuary-resistance  coefficient 

12.0 

PROPULSION  TESTS 

ref.  velocity  (knots)  RT  (lbs)  DeltaR 
23.7  ,  717.89,  191.1 

rpm  thrust  (lbs)  torque  (in-lbs) 

1175.0  ,  878.0  ,  2585.0 
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Percent  relative  uncertainties 

ref.  vel.  RT  DeltaR  rpm  thrust  torque 

1.4,  2.1,  2.1,  0.28,  0.71,  0.42 

Other  model  variables  :  slopes  of  1-t ,  JV  &  etaD  versus  CT 
dJV/dCT  d(l-t)/dCT  detaD/dCT 
-0.249  ,  0.067,  -0.015 

FULL-SCALE  VARIABLES  AND  UNCERTAINTIES 

Sea-water  kinematic  viscosity 
nu  X  10**5  (ft**2/sec) 

1.282 

Percent  relative  imcertainties  of  sea-water  properties 
density  kinematic  viscosity 
1.0,  2.0 

Full-scale  variables 

ship  length  (ft)  ship  speed  (knots) 

380.0  ,  25.0 

Percent  relative  uncertainties  of  full-scale  geometry 
length  area  prop  diameter 
0.5,  1.0,  0.1 

Percent  relative  uncertainties  of  full-scale  measurements 
speed  rpm  thrust  torque  SHP 
0.6,  0.4,  3.0,  0.9,  0.9 

SCALING  ALLOWANCE 
CA  dCA 
0.00035  ,  0.0001 
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EXAMPLE  OUTPUT  FILE 


LCCEXP 
DATE:  XXXX 
MODEL  No.  XXXX 
PROPELLER  No.  XXXX 
EHP  EXPERIMENT  No.  XXX 
SHP  EXPERIMENT  No.  XXX 
COMMENTS:  CASE  MFC 


INPUT  VARIABLES 


TANK-WATER  PROPERTIES 

water  density  (slug/ft**3) :  1.937 

kinematic  viscosity  coefficient  (ft**2/sec) :  0.1084E-04 
percent  uncertainty  of  density  :  0.14 

percent  imcertainty  of  viscosity  :  2.10 

MODEL  GEOMETRY 

length  of  model  (ft) :  22.697 

wetted  area  (ft**2) :  138.179 

diameter  of  propeller  (ft) :  0.9986 

percent  uncertainty  of  length :  0.14 

percent  imcertainty  of  wetted  surface  :  0.71 

percent  uncertainty  of  prop  diameter  :  0.07 

MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS 

residuary-resistance  coefficient :  0.00065 

percent  uncertainty  of  residuary-resistance  coef :  12.00 

MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS 

reference  velocity  (knots) :  23.70 

drag  RT  (lbs) :  717.89 

drag  DeltaR  (lbs) :  191.10 

propeller  rpm :  1175.00 

propeller  thrust  (lbs) :  878.00 

propeller  torque  (in-lbs) :  2585.00 

UNCERTAINTIES  OF  MODEL-SCALE  MEASUREMENTS 

percent  uncertainty  of  reference  velocity  (knots) :  1.40 
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percent  uncertainty  of  drag  RT  (lbs) :  2.10 

percent  uncertainty  of  drag  DeltaR  (lbs) :  2.10 

percent  uncertainty  of  propeller  rpm :  0.28 

percent  uncertainty  of  prop  thrust  (lbs) :  0.71 

percent  uncertainty  of  prop  torque  (in-lbs) :  0.42 


OTHER  MODEL-SCALE  VARIABLES 

slope  d  JV  /  d  CT :  -0.249 

slope  d  (1-t)  /  d  CT :  0.067 

slope  detaD/  dCT:  -0.015 

SEA-WATER  PROPERTIES 

kinematic  viscosity  coefficient  (ft**2/sec) :  0.1282E-04 
percent  uncertainty  of  density :  1.00 

percent  uncertainty  of  viscosity  :  2.00 

FULL-SCALE  SHIP 

ship  length  (ft) :  380.00 

ship  speed  (knots) :  25.0 

percent  xmcertainty  of  ship  length  :  0.50 

percent  xmcertainty  of  wetted  surface  :  1.00 

percent  xmcertainty  of  prop  diameter  :  0.10 

UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS 

percent  xmcertainty  of  ship  speed  :  0.60 

percent  xmcertainty  of  prop  rpm  :  0.40 

percent  uncertainty  of  prop  thrust :  3.00 

percent  xmcertainty  of  prop  torque  :  0.90 

percent  xmcertainty  of  SHP :  0.90 

SCALING  ALLOWANCE 

correlation  allowance :  0.00035 

xmcertainty  of  allowance  :  0.00010 

OUTPUT  VARIABLES 


MODEL-SCALE  VARIABLES 

total  resistance  coefficient  CT :  0.00246 

friction  resistance  coefficient  CF  :  0.00214 

residuary  resistance  coefficient  CR :  0.00065 
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advance  ratio  JV :  2.05 

thrust-deduction  factor  1-t :  0.60 

propulsive  efficiency  etaD  :  0.79 

UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES 

percent  uncertainty  of  CT :  4.14 

percent  uncertainty  of  CF :  0.37 

percent  uncertainty  of  CR  :  12.00 

percent  uncertainty  of  CA :  28.57 

percent  uncertainty  of  JV  :  1.43 

percent  uncertainty  of  1-t :  3.05 

percent  uncertainty  of  etaD  :  3.31 

FULL-SCALE  VARIABLES 

total  resistance  coefficient  CT :  0.00249 

friction  resistance  coefficient  CF  :  0.00149 

residuary  resistance  coefficient  CR  :  0.00065 

correlation  allowance  coefficient  CA  :  0.00035 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  SPEED 

percent  uncertainty  of  ship  speed  ;  0.60 

percent  uncertainty  of  prop  rpm :  1.60 

percent  uncertainty  of  prop  thrust :  6.90 

percent  uncertainty  of  prop  torque  :  6.57 

percent  uncertainty  of  SHP :  6.55 

percent  uncertainty  of  EHP :  5.57 

percent  overall  uncertainty  :  2.34 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  SHP 

percent  uncertainty  of  ship  speed  :  2.24 

percent  uncertainty  of  prop  rpm  :  2.62 

percent  uncertainty  of  prop  thrust :  5.16 

percent  uncertainty  of  prop  torque  :  2.79 

percent  uncertainty  of  SHP  :  0.90 

percent  uncertainty  of  EHP :  3.43 

percent  overall  uncertainty  :  1.37 
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APPENDIX  D: 

REPEATABILITY  OF  MODEL-SCALE  MEASUREMENTS 


The  precision  uncertainties  associated  with  submarine-model  large  cavitation 
channel  (LCC)  testing  are  investigated  by  considering  two  test  series.  These  two  test 
series  are  representative  of  the  submarine  model  resistance  and  powering 
experimental  evaluations  performed  in  the  LCC.  The  test  series  were  performed 
from  Jvme  through  September  1996.  The  current  model  test  speed  range  of  6-30 
knots  is  represented  in  these  tests,  as  well  as  the  current  data  collection 
instrumentation  and  calibration  techniques. 

The  precision  uncertainties  of  the  model  measurements  of  drag,  rpm,  thrust, 
and  torque  are  evaluated  in  two  ways.  First,  the  gage  calibrations  and 
instrumentation  were  analyzed  for  imcertainties.  The  drag,  thrust,  and  torque  gages 
are  calibrated  on  site  and  the  electronic  instrumentation  manufacturer  specifications 
are  examined.  These  precision  imcertainties  are  generally  very  small.  A  better 
assessment  of  the  model  measurement  precision  imcertainties  is  obtained  by 
analyzing  the  collected  test  data.  This  second  way  of  evaluating  the  measurement 
precision  uncertainties  takes  into  account  the  whole  data  collection  system, 
including  the  effects  of  changing  model  conditions  during  a  test  series,  water  flow 
variations,  variation  in  force  gages,  instrumentation  accuracy,  vibrations,  and 
computer  collection  and  recording  of  the  collected  model  drag,  shaft  RPM,  shaft 
thrust,  and  shaft  torque  values.  The  methods  used  to  determine  the  precision 
uncertainties  of  tiie  model  tests  provide  conservative  uncertainty  values  for  use  in 
the  global  uncertainty  analysis.  The  precision  imcertainties  for  the  measurements 
of  the  four  main  primary  model  quantities,  i.e.,  model  drag,  shaft  RPM,  thrust,  and 
torque,  are  now  examined. 

Calibration  of  the  drag,  thrust,  and  torque  gages  is  completed  before  each  test 
series.  The  instrumentation  currently  being  used  for  the  model  test  measurements 
is  presented  in  Table  D.l . 


A  typical  resistance  or  powering  experiment  consists  of  approximately  20  data 
spots  with  each  data  spot  representing  the  average  of  5  seconds  of  data  collected  at 
400  samples /second  for  a  total  sample  of  2000  for  a  set  speed,  drag,  and  RPM  (for 
powering)  condition. 

The  precision  uncertainty  in  model  drag  is  calculated  by  using  measured 
values  of  drag  at  the  same  nominal  model  speed.  A  correction  is  applied  to  the 
measured  drag  to  reduce  the  effects  of  the  speed  variation  on  the  measured  drag 
imcertainty.  The  measured  drag  is  multiplied  by  (V^)^  where  =  measured 

speed  and  V^,  =  nominal  speed.  EHP  tests  are  used  as  the  source  of  the  data.  The 
precision  imcertainties  are  presented  in  Tables  D.2  and  D.3  . 

The  precision  uncertainties  for  model  thrust,  torque,  and  RPM  are 
determined  differently  from  the  model  drag  uncertainty.  Typical  submarine 
powering  experiments  consist  in  varying  the  propeller  RPM  to  produce  different 
submarine  loadings.  The  precision  uncertainty  for  the  shaft  thrust,  torque  and 
RPM  has  been  estimated  by  determining  the  variation  of  the  thrust,  torque,  and 
RPM  data  from  a  least-square  curve  fit  through  the  data  spots  of  a  test.  Each  test 
contains  about  20  data  spots  which  comprise  a  range  of  thrust,  torque  and  RPM 
versus  total  drag  coefficient  (Q)  values.  The  thrust,  torque,  and  RPM  are  plotted 
against  and  a  second-order  least-square  curve  is  fitted  through  each  set  of  data. 

The  percent  difference  between  the  measured  data  spot  and  the  curve  at  each  Cj  is 
then  determined.  Twice  the  standard  error  estimate  (SEE)  divided  by  an  average 
thrust,  torque,  or  RPM  value  is  taken  as  the  precision  uncertainty  for  that  test.  The 
results  are  presented  in  Tables  D.4  and  D.5. 

The  following  table  summarizes  the  precision  uncertainties  for  LCC  model 
measurements  of  drag,  thrust,  torque,  and  RPM 


Uncertainties  of  model-scale  measurements  used  in  analysis 


drag 

thrust 

torque 

rpm 

1.5% 

0.5% 

0.3% 

0.2% 
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Table  D.1.  Instrumentation  used  at  the  LCC  for  resistance  and  powering  tests 
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ignal  Conditioner  Manufacturer  Valadyne  Scientific  Marine  Senries,  Inc 

Model _ CD19 _  IAF-01 


Table  D.2.  Model  drag  measurement  uncertainty  for  the  first 
LCC  resistance  and  powering  test  series 


1  0  knot 

16.878  ft/s 

1 8  knot 

30.380  ft/s 

RT  (ibs)** 
(data  spot) 

corrected  RT 
(10  knots) 

RT  (lbs)** 
(data  spot) 

corr  RT 
(18  knots) 

#  of  spots 

14 

14 

10 

10 

average 

102.539 

102.441 

308.413 

308.857 

2*Stdev 

1.289 

0.929 

5.109 

5.200 

% 

1.26 

0.91 

1.66 

1.68 

23.7  knot 

40.001  ft/s 

30  knot 

50.634  ft/s 

RT  (lbs)** 

corr  RT 

RT  (lbs)** 

corr  RT 

(data  spot) 

(23.7  knots) 

(data  spot) 

(30  knots) 

#  of  spots 

10 

10 

10 

10 

average 

520.428 

521.488 

823.579 

824.692 

2*Stdev 

2.514 

2.519 

18.680 

18.519 

% 

0.48 

0.48 

2.27 

2.25 

Total  Average  % 

1.33 

**  each  data  spot  =  an  average  of  2000  samples  collected  over  a 
5  second  collection  time  at  a  rate  of  400  samples/sec. 


Table  D.3.  Model  drag  measurement  uncertainty  for  the  second 
LCC  resistance  and  powering  test  series 


1 0  knot  1 8  knot 

16.878  ft/s _  _ 30.380  ft/s 


RT  (lbs)** 
(data  spot) 

corrected  RT 
(10  knots) 

RT  (lbs)** 
(data  spot) 

corr  RT 
(18  knots) 

#  of  spots 

12 

12 

12 

12 

average 

131.627 

132.370 

412.373 

413.243 

2*Stdev 

1.551 

1.560 

2.196 

2.187 

% 

1.18 

1.18 

0.53 

0.53 

23.7  knot 

40.001  ft/s 

30  knot 

50.634  ft/s 

RT  (lbs)** 
(data  spot) 

corr  RT 
(23.7  knots) 

RT  (lbs)** 
(data  spot) 

corr  RT 
(30  knots) 

#  of  spots 

12 

12 

12 

12 

average 

715.841 

717.771 

1169.741 

1170.641 

2*Stdev 

2.926 

2.658 

5.720 

4.614 

% 

0.41 

0.37 

0.49 

0.39 

Total  Average  %  0.62 


**  each  data  spot  =  an  average  of  2000  samples  collected  over  a 
5  second  collection  time  at  a  rate  of  400  samples/sec. 
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Table  D.4.  Model  thrust,  torque,  and  rpm  measurement  uncertainty  for  the  first 
LCC  test  series  at  model  test  speed  =  23.7  knots 
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Table  D.5.  Model  thrust,  torque,  and  rpm  measurement  uncertainty  for  the  second 
LCC  test  series  at  model  test  speed  =  23.7  knots 
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The  foregoing  table  defines  the  precision  uncertainties  for  model-scale 
measurements  of  drag,  thrust,  torque,  and  RPM  in  the  LCC.  The  uncertainties  in 
the  measurement  of  the  LCC  reference  flow  velocity  is  now  considered.  The 
reference  flow  velocity  for  model  tests  in  the  LCC  is  determined  by  means  of  LDV 
measurements  of  the  velocity  at  a  reference  point  in  the  flow  domain  where  the 
pressure  coefficient  vanishes.  The  location  of  this  reference  point  is  determined  via 
numerical  calculations.  Uncertainties  in  the  value  of  the  reference  velocity  stem 
primarily  from  two  main  sources  :  the  precision  errors  related  to  LDV  acquisition, 
and  errors  associated  with  the  location  of  the  reference  point  chosen  for  measuring 
the  reference  velocity. 

The  precision  error  associated  with  LDV  acquisition  is  the  dominant  cause  of 
uncertainty.  This  error  is  a  measure  of  the  turbulence  intensity  in  the  flow,  and 
thus  is  usually  minimal  for  free-stream  flows.  The  precision  error  stemming  from 
the  repeatability  of  LDV  acquisition  at  a  stationary  point  is  estimated  to  be  0.2%  . 
However,  this  error  increases  when  LDV  measurements  are  taken  at  different 
locations,  and/or  adjustments  are  made  to  the  optics.  Blanton  [3]  shows  that  the 
precision  error  associated  with  LDV  acquisition  can  be  as  high  as  0.8%  . 

The  secondary  contribution  to  the  uncertainties  of  the  LCC  reference  velocity 
stems  from  uncertainties  in  the  location  of  the  reference  point,  where  the  reference 
velocity  is  measured.  It  is  estimated  that  an  error  in  the  location  of  the  reference 
point  equal  to  one  inch  causes  a  0.1%  error  in  the  reference  velocity.  The  theoretical 
determination  of  the  reference  point  is  expected  to  provide  an  estimate  of  the 
location  of  the  reference  point  with  an  error  approximately  equal  to  1%  of  the  length 
of  the  model.  For  a  typical  20-foot  model,  the  location  of  the  reference  point  can 
then  be  presumed  to  be  known  with  an  error  of  2.5  inches.  In  addition,  LDV 
measurements  cannot  always  be  taken  exactly  at  the  theoretical  reference  points  due 
to  obstacles  and  restrictions  in  optical  access.  E.g.,  during  the  688  test,  measurements 
were  made  at  0.85  inch  off  the  theoretical  reference  point.  Obstacles  and  optical- 
access  difficulties  evidently  vary  from  model  to  model.  It  is  estimated  that  the 
distance  between  the  theoretically-determined  reference  point  and  the  reference 
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point  actually  used  in  the  LDV  measurements  can  be  as  large  as  4  inches.  The  total 
error  in  the  location  of  the  reference  point  where  the  reference  velocity  is  measured 
can  then  be  as  large  as  6.5  inches,  resulting  in  an  error  of  0.65%  as  was  explained 
previously. 

t  The  uncertainty  resulting  from  the  0.8%  precision  error  related  to  LDV 

acquisition  and  the  0.65%  error  associated  with  the  location  of  the  reference  point  is 

given  by  (  0.8^  +  0.65^  =  1.03%  .  Thus,  the  uncertainty  in  the  value  of  the 

reference  velocity  for  LCC  model  testing  is  taken  equal  to  1%  in  the  global 
uncertainty  analysis,  as  is  indicated  in  the  table  below. 

Uncertainty  of  LCC  reference  velocity  used  in  analysis 

reference  velocity 

1% 
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APPENDIX  E : 

UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS 

The  relative  total  (precision  +  bias)  uncertainties  of  full-scale  measurements 
of  the  ship  speed,  the  propeller  rpm,  the  thrust,  the  torque,  and  the  shaft 
horsepower  reported  for  four  full-scale  trials  (USS  Boise  SSN764,  USS  Colombus 
SSN762,  USS  Charlotte  SSN  766,  USS  Memphis  SSN  691)  are  given  below 


Reported  uncertainties  of  full-scale  measurements 


ship 

speed 

rpm 

thrust 

torque 

SHP 

SSN  764 

0.6% 

1.9% 

4.1% 

0.3% 

1.9% 

SSN  762 

0.6% 

0.5% 

4.1% 

0.3% 

0.6% 

SSN  766 

0.3% 

0.2% 

2.2% 

0.9% 

0.8% 

SSN  691 

0.2% 

2.2% 

1.4% 

0.9% 

Appreciable  variations  can  be  observed  in  the  foregoing  uncertainties.  Reasonable 
estimates  of  these  uncertainties  are  listed  below 


Uncertainties  of  full-scale  measurements  used  in  analysis 


speed 

rpm 

thrust 

torque 

SHP 

0.6% 

0.4% 

3.0% 

0.9% 

JO.9% 

These  estimates  of  full-scale  measurement  uncertainties  are  used  in  the  present 
uncertainty  analysis. 
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